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ABSTRACT 

Long gamma-ray bursts (GRBs) are associated with the explosion of massive stars in 
star forming regions. A large fraction of GRBs show intrinsic absorption as detected 
in optical spectra but absorption signatures are also detectable in afterglow X-ray 
spectra. We present here a comprehensive analysis the full sample of 93 GRBs with 
known redshift promptly observed by Swift XRT up to June 2009. The distribution 
of X-ray column densities clearly shows that GRBs are heavily absorbed indicating 
that they indeed occur in dense environments. Furthermore, there is a lack of heavily 
absorbed GRBs at low redshift (z 1 — 2) that might therefore be candidates for the 
missing 'dark' GRB population. However, there is no statistically significant correla- 
tion between the amount of X-ray absorption and the 'darkness' of a GRB. Finally 
we compare the hydrogen column densities derived in the optical with those derived 
from X-ray absorption. The two distributions are different, with the optical column 
densities being lower than the X-ray ones, which is even more apparent when correct- 
ing for metallicity effects. The most likely explanation is photoionization of hydrogen 
in the circumburst material caused by the radiation field of the burst. 
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1 INTRODUCTION 

(Long) gamma-Ray Bursts (GRBs) are recognized as being 
powerful cosmological tools, being detected from the local 
Universe to the highest redshifts reached today (z ~ 8.2) 
(Salvaterra et al. 2009; Tanvir et al. 2009). GRBs are so 
bright that they can illuminate, for a short time, regions of 
distant galaxies otherwise inaccessible and probe the cosmic 
chemical evolution of the Universe (e.g. Fiore et al. 2005; 
Prochaska, Chen & Bloom 2006; Fox et al. 2008). The GRB 
light shines through a medium that is very different (i.e. 
denser) from the typical interstellar medium (ISM) observed 
with traditional tools in distant galaxies High-resolution op- 
tical spectroscopic observations allow us to probe the rela- 
tive abundances of different elements highlighting the phys- 
ical state of the gas, its chemical enrichment and dust con- 
tent, and the nucleosynthesis processes in the GRB star- 
forming environment (Savaglio 2006; Prochaska et al. 2007). 
Despite the poor number statistics, a comparison between 
the metallicities obtained by GRBs and those from stud- 
ies of absorbing systems in the sightlines of quasar show 
that the former tend to be more metal-rich (on average by 
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a factor of 5 at high-,z) and to have a flatter or no red- 
shift evolution (see, e.g., Fynbo et al. 2008, 2009; Savaglio, 
Glazebrook & LeBorgne 2009; Chen et al. 2009). Variations 
in the intensity of optical absorption lines were detected in 
a few cases of very optically bright GRBs. These variations 
help to constrain the distance of the absorbing gas to the 
GRB. This has been done in only two cases where distances 
larger than ~ 2 kpc have been derived. Matter within this 
radius has been photoionized by the GRB radiation flux, 
which means that optical observations can only probe the 
outermost part of the region surrounding the GRB (D'Elia 
et al. 2009; Vreeswijk et al. 2007). 

Matter along the line of sight also attenuates X-ray 
photons. X-rays are photoelectrically absorbed by metals, 
producing a well known and well identifiable pattern in the 
observed X-ray spectra (Morrison & McCammon 1983). In 
contrast to absorption in the optical range, the absorption 
of X-rays is not sensitive to the state of the element (either 
gas or solid) , therefore these measurements provide an unbi- 
ased view of the total absorbing column density of material. 
As demonstrated by several studies the X-ray absorbing col- 
umn towards GRBs is high (Stratta et al. 2004; Campana et 
al. 2006; Watson et al. 2007), which implies in principle, i.e. 
for any known dust to gas ratio, a high dust extinction. In 
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contrast, GRB optical afterglows, when present, often show 
a lack of reddening. Several suggestions have been advanced 
to explain this dichotomy but a clear understanding has not 
been reached yet (e.g. Galama & Wijers 2001, see also Nar- 
dini et al. 2009). 

Stratta et al. (2004) presented a systematic analysis of a 
sample of 13 bright afterglows observed with the BeppoSAX 
narrow field instruments. In only two cases they found a sig- 
nificant detection of additional intervening material in ex- 
cess of the one in our Galaxy. However, due to the limited 
photon statistics, they could not exclude that intrinsic X-ray 
absorption was also present in the other bursts. Campana 
et al. (2006) presented a systematic study of the intrinsic 
column density absorption in 6 GRBs with known redshift 
observed by Swift. Adding these GRBs to others found in the 
literature with known redshift (mainly Chandra and XMM- 
Newton observations), they concluded that there is strong 
evidence that GRBs are intrinsically absorbed (Campana 
et al. 2006). The distribution of absorbing column densities 
is consistent with GRBs occurring within molecular clouds, 
which supports the link between (long) GRBs and star for- 
mation. Watson et al. (2007) found a large discrepancy be- 
tween optical and X-ray derived GRB column densities and 
suggested photoionization of the gas cloud as the most likely 
cause. 

Here we consider the full sample of 93 GRBs detected 
by Swift (Gehrels et al. 2004) XRT (Burrows et al. 2005) 
up to June 2009 which were observed within 1,000 s (to 
ensure relatively high signal) from the burst onset and have 
a known redshift. In section 2 we briefly describe the data 
analysis and the way we evaluated the column density. In 
section 3 we discuss our results and in section 4 we draw 
our conclusions. 



2 DATA ANALYSIS 

We made extensive use of automated data products pro- 
vided by the Swift/XRT GRB spectra repository (Evans et 
al. 2009). For each burst all the data are processed in an 
homogeneous way and summed together. Most GRBs have 
one spectrum for each of the two instrument modes (Win- 
dowed Timing, WT, and Photon Counting, PC; Hill et al. 
2004) . For each time interval, a source and background spec- 
trum, an exposure map and an ancillary response file (arf) 
are created and then combined together. To identify inter- 
vals affected by pile up we search for time intervals where 
the count rate is above 0.6 counts s _1 in PC mode or 150 
counts s _1 in WT mode. For those intervals we then obtain 
the profile of the source and compared it to the calibrated, 
non-piled-up, Point Spread Function (Moretti et al. 2005), 
for both PC and WT mode separately. Background spectra 
are created in an analogous way excluding contaminating 
sources in the background region. The individual arf files, 
which are generated in order to account for CCD defects, 
are combined using the addarf tool, weighted by the counts 
in the individual source spectra. 

For each mode one source and one background spec- 
trum are created. The source spectrum is binned using the 
grppha tool such that there is at least one count in each 
spectral bin necessary for XSPEC to correctly calculate the 
Cash-statistic (K. Arnaud, priv. communication) used to fit 
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Figure 1. Distribution of the X— ray column densities of Swift 
GRBs promptly observed by the XRT instrument. The continu- 
ous line histogram has been derived from positive detections, the 
dashed histogram also includes upper limits. The continuous line 
represents the fit with a lognormal function. 



the spectra (Cash 1979). The spectra are modelled with 
an absorbed power law within XSPEC12. The absorption 
model adopted is phabs. The abundances are fixed at the 
solar value from Anders & Grevesse (1989). The photoelec- 
tric absorption consists of two components: one from our 
Galaxy, kept constant to the value taken from Kalberla et 
al. (2005), and the other one free to vary, shifted in energy 
at the (fixed) redshift of the GRB. 

Evans et al. (2009) provided already a collection of col- 
umn densities of Swift GRBs. However, either they did not 
fix the redshift in the fit of the X-ray column density and, 
when they fix it, they also included bursts without known 
redshift fixed at the mean value of the sample. In addition, it 
has been shown that the absorbing column density can vary 
significantly when fitting an afterglow X-ray spectrum with 
a power law model when the overall spectrum is (slowly) 
changing (Butler & Kocevski 2007; Falcone et al. 2006). 
In the case of a slowly varying spectrum a cut-off spectral 
model can be a solution to this problem. 

For our sample, we selected GRBs with known spectro- 
scopic redshift and with observations starting within 1,000 s 
from the burst onset. The knowledge of the redshift is nec- 
essary in order to derive the intrinsic column density. With 
these criteria our sample contains 93 GRBs up to May 2009. 
For these GRBs we derived the intrinsic absorbing column 
density selecting the time interval in the light curve when 
the 0.3-1.5 keV to 1.5-10 keV hardness ratio is constant. 
During the first fast steep decay time interval the spectrum 
is strongly variable, as observed in several bursts. Therefore 
this interval is left out as well as any interval with (major) 
flaring activity. With these criteria we preferentially select 
observations in the 3, 000 — 10 5 s time interval and in PC 
mode, after the early steep decay phase (Nousek et al. 2006). 
In a few cases when the late-time signal was not high enough 
to derive a constraint on the intrinsic column density, we se- 
lected short intervals in the early light curve during which 
the hardness ratio was constant (usually in WT mode). A 
comprehensive report of the selected time intervals and re- 
sults is provided in Table 1. 
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Table 1. Column densities for Swift GRBs with known rcdshift. 
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Table 1. continued. 
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Errors at 90% confidence level. Upper limits are at 3cr level. For the error on the power law slope we just quote the maximum between 

lower and upper values. 

PC and WT refer to the Swift XRT observing mode. SM or VAR is a small or variable hardness ratio. D3 and B3 refer to burst which 
are optically dark or bright at 1,000 s. Data arc taken from Cenko et al. (2009). D3 are bursts with the optical (R) to X-ray spectral 
index /Sox evaluated at 1,000 s lower than 0.5 (Jakobsson et al. 2004). B3 bursts have f3ox > 0.5. Dll and Bll are dark and bright 
bursts from the sample of Zheng et al. (2009) where Pox has been evaluated at 11 hr. Bll are bursts with we took fiox > 0.6. Note 

that some GRBs that are 'dark' at 1,000 s are considered 'bright' at 11 hr. 
The last three GRBs in the table have XRT observations starting later than 1,000 s. We list them here because they have an optically 

determined hydrogen column density. 



3 RESULTS AND DISCUSSION 

We analyzed the XRT spectra of 93 GRBs promptly ob- 
served by Swift with known redshift (up to May 2009). For 
only 8 of them we do not find evidence of intrinsic X-ray 
absorption at the host galaxy site. This clearly indicates 
that the medium surrounding GRBs is not tenuous. In con- 
trast to optical studies which are subject to uncertainties 



due to ionizations level and elements depleted to dust, X- 
ray absorption studies are less biased since the absorption is 
produced by the innermost electrons and therefore does not 
depend on the state of the absorbing medium. The distribu- 
tion of intrinsic absorbing column densities is shown in Fig. 
1. This distribution can be well described by a lognormal 
distribution with mean \ogN H {z) = 21.9 ± 0.1 (90% confi- 
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Figure 2. X-ray column densities versus redshift. Upper limits 
also plotted as arrows. 



dence level) and standard deviation <J\ og N H {z) = 0.5 ± 0.1 
(the fit is good with a reduced % 2 of 1.2). Given the size of 
our sample this clearly indicates that X-ray afterglow spec- 
tra are intrinsically absorbed (see also Nardini et al. 2009). 
The amount of this absorption is typical of the plane of our 
Galaxy considering solar metallicity. 

In order to determine the origin of the absorption ex- 
cess, we compare the distribution of measured intrinsic col- 
umn densities with the distribution expected for bursts 
occurring in different ambient media. We already demon- 
strated (Campana et al. 2006) that the column density dis- 
tribution is consistent with Galactic-like molecular clouds 
(Reichart & Price 2002) and not with a random GRB oc- 
currence within a spiral galaxy like our own (Vergani et 
al. 2004). Our extended sample strengthens this conclusion 
since a random occurrence in our Galaxy would predict a 
sizable (~ 30%) GRB population with no intrinsic absorp- 
tion which is not present in our sample. Following Reichart 
& Price (2002) we parametrize the column density distribu- 
tion with a broken power-law 



n(log Nh) 



he (N H ) b 



x-rn-n— )" (0<N H <a) 
) C (Na>a) 



1 he I Nh 
In 10 c-b \ a 



(1) 



where a is the break column density, b is the low column 
density power-law index, and c is the high column density 
power-law index. We find that log (a/cm -2 ) = 21.7li !i5) 
b = 1. 59^0 57, and c = — 0.78±g;^ (90% confidence level). 
We stress that our sample also contains upper limits and 
that we are sensitive to low values of column density, which 
however are found only in a small fraction of the total sam- 
ple. If we include 6 additional tight upper limits (limits with 
Nh < 4 x 10 21 cm -2 ) we obtain very similar values as ex- 
pected because these additional GRBs comprise just 7% of 
the sample. These results strongly support an origin of long 
GRBs within high density regions of host galaxies. 

With our large sample we can also check if there is 
any dependence on redshift. To detect an intrinsic column 
density, one preferentially selects higher column density at 
higher redshifts since the column density contribution is 
shifted towards lower energy bands and therefore more dim- 



Figure 3. X-ray column densities versus optically derived col- 
umn densities along with their errors. Dashed lines indicate val- 
ues within a factor of two from the line of equal X-ray to column 
density (upper continuous line). Dotted lines mark optical column 
densities nth orders of magnitude less than X-ray ones. 



cult to identify. To get the same absorption contribution in 
an X-ray spectrum at different redshift, the intrinsic column 
density should increase as ~ (1 + z) 2,6-2 ' 7 (e.g. Galama & 
Wijers 2001). The interesting point is that of the ten GRBs 
at redshift z £ 4 only two have upper limits (see Fig. 2). 
The mean value of this absorption (excluding upper limits) 
is ~ 5 x 10 22 cm" 2 with a large dispersion. This indicates 
that the available space for a weakly absorbed GRBs is at a 
level of only ~ 20%. To quantify this statement we compare 
the column density distribution of bursts at z > 4 with the 
ones at z < 1. For the two upper limits in the z > 4 sam- 
ple we take a worst case approach an fix them to the mean 
column density of the z < 1 sample (i.e. 5 x 10 21 cm -2 , see 
also below). A Kolmogorov-Smirnov (KS) test of the two 
distributions shows that they are not drawn from the same 
parent distribution with a probability of 0.08% (equivalent 
to Gaussian 3.3 <r). This value does not change for lower 
values of the high-z upper limits. This implies either the 
lack of low absorption burst at high redshift or the lack of 
heavily absorbed bursts at low redshift. The first instance 
might be explained by a higher mean and/or compactness 
of star forming regions at high redshifts. Alternatively, in 
a few cases, absorption caused by intervening systems can 
also play a role (Campana et al. 2006). The second instance 
might find an explanation in the change of optical extinc- 
tion curve (e.g. due to a change in dust composition, dust to 
gas ratio, grain size): at high redshifts the extinction curve 
should be much flatter and similar to the one determined for 
quasars (Maiolino et al. 2004), at lower redshifts the curve 
should resamble more closely the one of our Galaxy, with a 
higher optical absorption for the same X-ray absorbing col- 
umn density. If this is true heavily absorbed bursts at low 
redshift are missing from our sample because for these bursts 
we cannot get an optical counterpart and/or a redshift, i.e. 
they likely populate the 'dark' burst class. 

We therefore search for a possible correlation between 
the presence of high X-ray absorbing column density and 
the darkness of a GRB. Following Jakobsson et al. (2004) 
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a GRB is classified as dark if the optical (R band) to X- 
ray spectral index Pox is smaller than 0.5. We first consider 
the sample of Cenko et al. (2009) where Pox has been evalu- 
ated at 1,000 s after the burst detection (see also Perley et al. 
2009). In this sample there are 5 dark bursts and 9 non-dark 
bursts common to our sample. A KS test on the distribution 
of column densities of the two sample provides a probability 
of 15% of being drawn from the same population. We also 
consider the work by Zheng et al. (2009) in which Pox has 
been evaluated at 11 hr (as the original definition by Jakobs- 
son et al. 2004). The two optical samples are not overlapping 
and some burst classified as 'dark' in one sample is not dark 
in the other and viceversa. In this case we have 4 dark GRBs 
and 39 firm non-dark GRBs (with Pox > 0.6) in common 
with our sample. We included in the non-dark burst sample 
three upper limits (with Nh < 4 x 10 21 cm~ 2 ) taking their 
upper limit as the true value. In this case a KS test provides 
a probability of 2% (corresponding to a 2.3 a Gaussian prob- 
ability). This might suggest a different population but the 
dark population lack a sufficiently high number of GRBs to 
perform a reasonable statistical analysis. 

It has been suggested that not only intrinsic absorption 
might be a cause of darkness of a GRB but also the redshift. 
We therefore tried a two parameters (or bi-dimensional) KS 
test (Press, Flannery & Teukolsky 1986), testing together 
the X-ray column density and the redshift. In the case of 
Pox evaluated at 1,000 s, the 2D-KS probability is 16%, 
in the case of Pox evaluated at 11 hr the 2D-KS probabil- 
ity is 16%. This indicates that X-ray absorption alone or 
together with redshift cannot provide a statistically sound 
explanation for the darkness of GRBs. Cenko et al. (2009) 
suggested that optical absorption is the likely cause, but 
based on their sample a KS test of the host absorption of 
dark and non-dark GRBs cannot statistically support this 
suggestion. 

Last, we compare X-ray column densities to the hydro- 
gen column densities derived from optical data (available 
for GRBs at z ^ 2). This is now possible thanks to our large 
X-ray sample and the optical GRB sample recently pre- 
sented by Fynbo et al. (2009). These two samples comprise 
36 common GRBs (including three more from the literature, 
GRB050505, GRB050904, and GRB081203). Out of these 36 
GRBs, 5 GRBs have only an upper limit in X-rays and one 
has an upper limit in the optical (see Fig. 3). A KS test of 
the 30 GRB X-ray and optical column densities provide a 
probability of 1.1% (Gaussian 2.5 a) that the two popula- 
tions come from the same distribution. As we see in Fig. 3, 
the bulk of the GRBs have column densities in the X-rays 
which are within a factor of ~ 10 larger than the ones in 
the optical. In addition, none of the bursts has an optical 
column density larger than a factor of ~ 2 compared to the 
X-ray one, including errors, which indicates that X-ray col- 
umn densities are, as expected, good tracers of the ongoing 
absorption. This is not the case for optical column densities 
since a large number of them lies well below the X-ray de- 
rived value. The most extreme GRBs have ratio of X-ray 
to optical column density in the 10~ 5 — 10~ 4 range. Several 
possibilities have been discussed in the literature. For differ- 
ences of an order of magnitude, higher metallicities might 
be invoked, providing more metals (X-ray absorption) at 
the same value of hydrogen (Campana et al. 2007, 2008; 
Watson et al. 2007). For larger differences the likely expla- 
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Figure 4. X-ray versus optical column densities. X-ray column 
densities were corrected for metallicities different from the solar 
ones adopting the ones measured in the optical. The continuous 
line refers to equal optical to X-ray values, the dotted lines to 
nth orders of magnitude difference. 

nation relies on ionization of hydrogen by the GRB high 
energy flux (Lazzati, Perna & Ghisellini 2001; Watson et al. 
2007) in the close GRB environment. 

From spectroscopic studies of quasar Damped Lyman 
Alpha (DLA) we know that the Sill A1526 line intensity can 
be used as a proxy for metallicity M through the relation 
[M/H] = -0.92 + 1.41 Iog(£W/A) (with EW being the 
equivalent width of the Sill line; Prochaska et al. 2008). 
GRBs seem to follow the same correlation, although possibly 
with a somewhat steeper slope, but the sample of bursts 
with optically derived metallicities and the detection of the 
Si II A1526 line is still very limited. Following this suggestion 
we derived a rough metallicity value for 29 GRBs based on 
the Sill line strength. For 24 of them (including 3 upper 
limits) we have also a determination of the optical column 
density. With this sample we first investigate if there is a 
correlation between the derived metallicity and the X-ray 
column densities and find that there is none. As a further 
step we assume that the metallicity derived from the Sill 
lines applies also to the medium absorbing the X-ray flux. 
If this is the case the true column density absorbing the X- 
ray flux is just the derived value divided by the metallicity 
value, since the X-ray fitting only allows to determine the 
integrated contribution of Nh /M. We can compare in Fig. 
4 the 'corrected' X-ray column densities with the optically 
derived ones for 24 GRBs. In this case no metallicity effect 
should be present. The disagreement between the two is even 
more severe. The probability that the distribution of the 
optical and X-ray column densities were drawn from the 
same parent population is 8 x 10 -7 , equivalent to Gaussian 
4.9 a, which calls for ionization of hydrogen by the GRB 
high energy flux. 



4 CONCLUSIONS 

We considered a complete sample of 93 Swift long GRBs 
with prompt XRT pointing ( ^ 1,000 s) and with a spec- 
troscopic redshift. The distribution of the intrinsic X-ray 
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absorption column density is consistent with a lognormal 
distribution with mean log N H (z) = 21.9 ± 0.1 (90% confi- 
dence level) and standard deviation cr\ og N H (z) = 0.5 ± 0.1. 
This confirms that long GRBs occur in a dense environment, 
which is consistent with the expected distribution of absorp- 
tions for GRBs occurring randomly in molecular clouds (and 
not with long GRBs occurring randomly in any place of a 
galaxy similar to own Galaxy). 

Looking at the distribution of X-ray column densities 
versus redshift we find a lack of non-absorbed GRBs at high 
redshift and a lack of heavily absorbed GRBs at low red- 
shift. The former might be explained in terms of more com- 
pact and dense star formation regions in the young Universe 
(or to a sizable contribution from intervening systems) . The 
latter might be interpreted as due to a change in the dust 
properties with redshift, with GRBs at redshift z £ 2 — 3 
having a higher dust to gas ratio for the same X-ray col- 
umn density. This will naturally provide a lack of heavily 
(X-ray) absorbed GRBs at small redshifts, i.e. dark GRBs. 
We also investigated in more detail the relation between X- 
ray absorption and darkness of a GRB. We considered two 
different samples of long GRBs selected at different times 
(1,000 s and 11 hr). For none of the two we found a sta- 
tistically significant relation between X-ray absorption and 
darkness, nor in a two-dimensional test including X-ray ab- 
sorption and redshift. We note however that the sample of 
dark GRBs with known redshift is small and the weakness 
of our correlation might also depend on this. 

Finally, we compare the hydrogen column densities de- 
rived in the optical and those from X-ray absorption. The 
two distributions are markedly different, with the optical col- 
umn densities being lower than X-ray ones (and even more 
if correcting for metallicity effects). The most likely expla- 
nation is in terms of photoionization (depletion) of hydrogen 
in the circumburst material caused by the burst radiation 
field (Lazzati et al. 2001; Watson et al. 2007). 
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